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ABSTRACT 

We present an X-ray eclipse timing analysis of the transient low mass X-ray bi- 
nary XTE J1710-281. We report observations of 57 complete X-ray eclipses, spread 
over more than a decade of observations, made with the proportional counter array 
detectors aboard the RXTE satellite. Using the eclipse timing technique, we have 
derived a constant orbital period of 0.1367109674 (3) d, during the period from MJD 
52132 upto MJD 54410; and la limits of -1.6 x 10 j ^-i and 0.2 x 10 j 
d~^, on the period derivative, Porb- This puts constraints on the minimum timescale 
of secular orbital period evolution {Porb/ Porb) of 2.34 x 10® yr for a period decay 
and 18.7 x 10® yr for a period increase. We also report detection of two instances of 
discontinuity in the mid-eclipse time, one before and one after the above MJD range. 
These results are interpreted as three distinct epochs of orbital period in XTE J1710- 
281. We have put lower limits of 1.4 ms and 0.9 ms on orbital period change (APorh) 
at successive epochs. The detection significance of the two orbital period glitches are 
llfT and 4(7 respectively. The sudden changes in orbital period is very similar in nature 
to that observed in EXO 0748-676, though their magnitude is much smaller in XTE 
J1710-281. 

Key words: accretion, accretion discs, binaries: eclipsing, binaries: general, stars: 
individual: XTE J1710-281, stars: neutron. X-rays: stars 



1 INTRODUCTION 

XTE J1710-281 is a transient Low Mass X-ray Bi- 
nary (LMXB) which was discovered in 1998 by the 
Rossi X-ray Timing Explorer [RXTE), and is likely to 
be assoc iated with the ROSAT source IRXS J171012.3- 
280754 (|Markwardt et al.l Il998l '). It is a highly variable 
source and several bu r sts have been repor ted to occur 
iMarkwardt et all I2001I : iGallowav 6^111120081'). The system 
has an orbital period of 3.28 hr ijMarkwardt et al.l I2OOII ) 
and the hght curve shows dipping phenomena which could 
be due to occultations in the outer regions of the accre- 
tion disk as seen in m any other high inclination LMXBs 
(|White fc Swanklll982l ). 

XTE J1710-281 is a poorly studied LMXB. Being 
in a binary system, the orbital period of XTE J1710- 
281 is expected to change, due to redistribution of the 
angular momentum arising from interaction between the 
components of the binary system. The orbit can evolve 
due to various mechanisms, such as, mass transfer within 
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the system due to Roche lobe overflow, tidal interaction 
between the components of the binary system, gravita- 
tional wave radiati on, magnetic braking, (|Rappaport et akl 
19831: Hurley et al. 2002 ) and X-ray irradiated wind outflow 
( Ruderman et al.ll 19891 ). The measurement of orbital period 
derivative and hence the orbital evolution is therefore, im- 
portant to understand the physical processes occuring in the 
system. 

Orbital evolution can be determined by several tech- 
niques. Measurement of pulse arrival time delay, is one of the 
methods used to de t ermine the evolu t ion of the binary orbits 
(|Deeter et al.lll99ll : [Paul et all bool . 120071 ). The method of 
pulse folding and x'^ maximization with a varying orbital 
ephem eris, has also been used to determ ine the orbital evo- 
lution (|Paul et al.ll2002l : I Jain et al.ll2007l '). Both these meth- 
ods are well established and statistically at par. However, 
to determine the pulse arrival times, the data-length for 
each sample should be kept small so as to avoid significant 
smearing of pulse phase due to orbital motion. In the case 
of faint sources, or in the case of observation made with a 
small photon collection area (e.g. with focussing optics), a 
small integration time may be insufficient for measurement 
of pulse arrival time. In such a case, maximisation of pulse 
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detection by varying the orbital ephemeris is an effective 
techn ique as the entire data set is used together l|Paul et alj 
Eclipse timing is another technique, u sed to deter- 
mine the orbita l evolution of LMXBs ( Wolff et all l2009l : 
I Jain et aljboiol ). In some sources, in the absence of pul- 
sations or eclipses, orbital period derivative has been mea- 
sured with s ome stable orbital intensi ty modula tion features 
f4U 18 20-30: IChou fc Grindlavl l|200ll ): Cyg X-3: ISingh et~al] 
I 2002l )l. 

XTE Jf7f0-28f is one of the very few LMXBs, 
where full, sharp X-ray eclipses have bee n observed. 
The other systems b eing, EXO 0748-676 (|Wolff et all 
120021 ) GRS J1747-312 Jin't Zand et al.ll2003l ). MXB 1658- 
298 jCominskv fc WoodI Il989l ) and AX J 1745. 6-2901 
l|Maeda et all Il996l). Among these sourc es, EXO 0748- 
676 (|Parmar et all 1 19861 : 1 Wolff et al.ll2002l ) is the only sys- 
tem in which a large number of eclipse s have been timed 
with high accuracy (|Wolff et all \200^ . In case of GRS 
J1747-312 and AX J1745.6-2901, the eclipse duration is 
too long (^43 minu t es and ^23 minutes , respectively, 
lin't Zand etal] (|2003t ): iPorguet et ai] l|2007l ')) to carry out 
monitoring measurements. AX J1745. 6-2901 is located near 
a bright source, which puts strong constraints on the timing 
analysis. The fifth source, MXB 1658 -298 has been mostly 
in inactive state, since its discovery l|Wachter et al.l l2000l : 
lOosterbroek et al.ll200lh . The LMXB XTE J1710-281 has a 
fairly small eclipse duration and has been persistently active 
since its discovery. This makes it an ideal source to investi- 
gate the orbital evolution. 

In this work, we have used the sharp eclipses of XTE 
J1710-281 as timing markers to determine the orbital period 
of the system. By measuring the mid-eclipse times, over a 
long time baseline of ~11 years, we can determine the change 
in the orbital period and hence estimate the orbital evolution 
of the system. 



2 OBSERVATIONS AND ANALYSIS 

Data for the present analysis were obtained from ob- 
servations made with the Proportional Counter Array 
(PGA) oi l board the Rossi X -ray Timing Explorer (RXTE) 
satellite jBradt et al.l Il993l 'l. The RXTE-PCA consists of 
an array of five coUimated xenon/methane multianode 
proportional counter units (P GU) with a total photo n 
collection area of 6500 cm^ (|jahoda et all 1 19961 . l2006h . 
depending on the number of PGUs ON. The entire analysis 
was done using ftools from the astronomy software 
package HEASOFT-ver 6.10. The PGA data collected in 
the event mode and the Good Xenon mode, were used 
to generate the light curves, using the ftool-SEEXTRCT. 
The analysis was done in the energy band 2—20 keV. The 
background was estimated using the ftool-pcabackest. 
Faint source model was taken from the RXTE website 
(http:/ /heasarc. gsfc.nasa.gov/docs/xte/pca_ news.html). 
Thereafter, the background subtracted light curves were 
barycenter corrected using the ftool-faxbary. 

We have analyzed all the RXTE-PCA archival data, 
which covered a full X-ray eclipse. However, since this is a 
bursting source, we ignored few eclipses, where bursts oc- 
curred close to the ingress and egress phase of the eclipse. 
From data spread over ~11 years (1999-2010), we have found 
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Figure 1. A sample of 2—20 keV, background subtracted light 
curve of XTE J1710-281, binned with 2 s and including an eclipse 
lasting for about 428 s. The line represents the best fit five- 
parameter model for the light curve. The best fit had a of 
381 for 375 d.o.f. 



57 complete eclipses. Table 1 gives the observation IDs of all 
the 57 eclipses, observed with RXTE-PCA. In Figure 1, we 
have shown a sample background subtracted light curve of 
XTE J1710-281 (Obs ID: 91045-01-01-13), binned with 2 
seconds and including an eclipse lasting for ~ 420 s, exclud- 
ing the ingress and egress phase. 

For most of the observations, apart from a few type-I 
X-ray bursts, the out-of-eclipse count rate did not seem to 
have any significant variability. Therefore, the variable com- 
ponents of the light curve around the eclipse phase are, the 
pre-ingress (Cpre-ingress), eclipse (CecUpse), and post-egress 
(Cpost-egress) couut ratcs; the ingress (Ati„) and egress 
(Ateg) duration, the eclipse duration (Ateci) and the mid- 
eclipse time. Considering all the components to be freely 
variable, we first fitted a seven-paramete r ramp a nd step 
model to the eclipse phase (similar to Wolff et al.l |2009j)). 

pre — ingress and Cpost — egress 



It was found that the value of C 



were similar and the eclipse ingress and egress duration were 
also similar withing errors. The parameter space, thus got 
reduced to five: the pre-ingress and the post-egress count 
rate (Cpre-in-post-eg), the eclipse count rate, the ingress 
and egress duration {Ati„^g), the eclipse duration (Ateci), 
and the mid-eclipse time. The average eclipse duration in 
XTE J1710-281 is ~ 420 s; and for the model fitting, ~ 150 
seconds of data was taken before and after the eclipse phase. 
It was also seen that the error in the mid eclispe time mea- 
surement was smaller when the five-parameter model was 
used. The best fit model for the eclipse profile in Figure 1, 
is shown with a solid line. The best fit had a of 381 for 
375 degrees of freedom. 

All the 57 X-ray eclipse light curves were fitted with a 
five-parameter ramp function, as described above. The mid- 
eclipse times and the corresponding la errors were deter- 
mined. The results are given in Table 1. The errors in the 
mid-eclipse times vary between 0.000006 d to 0.000064 d, i.e. 
0.5 s to 5.5 s. This is mainly due to difference in the relative 
count rates of the source and the background, and partly 
due to the number of detectors ON. The orbit numbers are 
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Figure 2. The observed minus calculated times (residuals, O- 
C) for eclipses observed in XTE J1710-281 during 1999-2010, 
obtained from RXTE-PCA observations. The O-C variation is 
plotted as a function of orbital cycle. The three distinct orbital 
period epochs are mentioned. 



with respect to the first echpse detected in the RXTE-PCA 
data. We fitted a constant and a linear model to the eclipse 
measurements between MJD 52132 - 54410. The results of 
the fit are given in Table 2. We obtained an orbital period of 
0.1367109674 (3) d (epoch MJD 51250.924540 (4)) and have 
derived la limits of 0.2 x 10 "^^ d d"^ and -1.6 x 10 "^^ d 
d^^, on the orbital period derivative (Port). The of the 
fit was 74 for 51 d.o.f. Before and after the above mentioned 
MJD range, we found shifts in the mid-eclipse times. We 
refer to these shifts as three epochs in the orbital period. 

Figure 2 shows the "observed minus calculated" (O - C) 
diagram for all the eclipse measurements of XTE J1710-281, 
obtained after subtracting the linear component obtained 
from epoch 2. The three different epochs of orbital period 
is evident. It is obvious from the figure that a polynomial 
function consisting of linear (Port), quadratic (Porb), cubic 
{Porb) etc terms cannot be fitted to the observed dataset. 
A piecewise linear function could be more appropriate. But, 
there are few observations in epoch 1 and 3, hence one can- 
not determine the orbital period during epoch- 1 and epoch-3 
with very high accuracy. However, from our observations, we 
have put lower limits on orbital period changes of AP — lA 
ms (1.7 X 10 ~* d) between epoch-1 and epoch-2; and a AP 
= 0.9 ms (1.1 X 10 ~* d) between epoch-2 and epoch-3. The 
detection significance of the two orbital period glitches are 
11a and ia respectively. 

We have also created average eclipse profiles for the 
three epochs by combining all the corresponding eclipse light 
curves. The eclipse light curves were co-added using the or- 
bital period and epoch given in Table 2. The folded pro- 
files are shown in the top panel of Figure 3. The best fit 
ramp function is also shown for the eclipse profile of epoch 
2. The best fit had a of 689 for 650 de grees of freedom 
(d.o.f). The bottom panels of the same figure shows the 
eclipse ingress and egress profiles. It is clear from Figure 
3 that the time of ingress and egress of eclipse, during the 



three epochs is shifted in phase while the eclipse duration 
has remained nearly same. 



3 DISCUSSION 

Measurements of change in orbital period of an LMXB, is 
a crucial diagnostic to understand the accretion processes 
occuring in a binary system; and their effect on the system 
parameters. And the orbital period can be very well deter- 
mined by time connecting a stable fiducial marker in the 
light curve of the binary system. We have analyzed 57 full 
X-ray eclipses of XTE J1710-281, observed by the RXTE 
satellite. The observations cover more than 30000 binary or- 
bits spread over ~11 years. A five-component model was 
fit to each eclipse profile and the mid-eclipse times and the 
corresponding errors determined. We have determined an 
orbital period of 0.1367109674 (3) d and limits on the pe- 
riod derivative of -1.6 x 10 "^^ d d"^ and 0.2 x 10 "^^ d 
d~^, i.e., limits on the timescale of secular orbital period 
evolution, Porb/ Porb, of 2.34 X 10 * yr for a period decay 
and 18.7 x 10 * yr for a period increase, respectively, during 
the period from MJD 52132 to MJD 54410. 

The variation in the orbital ephemerides of 
XTE J1710-281, is significantly different from that 
seen in m ost of the other LMXB systems, such as, 4U 
1820-303 jChou fc Grindlayl 1200 ll). SAX J1808.4-3658 
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20071). Her X -1 (|Paul et alll2007l). X 2127-H1 9 



Homer fc Charleslll998l ) and 4U 1822-37 (|jain et al.ll201(il ') 



During the period from MJD 52132 to MJD 54410, the 
limits on the orbital period derivative of XTE J1710-281, 
is more than an order of magnitude smaller than those 
measured in the other LMXBs (4U 1820-303 (3.5 x 10 
~* yr"^), SAX J1808.4-3658 (1.3 x 10 yr'^). Her X-1 
(-2 x 10 yr-i), X 2127-^119 (9 x 10 yr"^) and 
4U 1822-37 (2.0 x 10 "'^ yr"^)). Outside the above MJD 
range, the observed trend in the residual (O - C) behaviour 
of XTE J1710-281, is also different from that seen in the 
aforementioned LMXBs. 

The observed O - C variation strongly resemble the one 
seen in EXO 0748-676 l| Wolff et al.ll2009i ). Interestingly, of 
the known LMXBS with well determined eclipse times, EXO 
0748-676 and XTE J1710-281 have shortest duration of the 
eclipse, making it easier to monitor with X-ray observa- 
tories. Though fewer eclipses have been observed in XTE 
J1710-281, as opposed to more than 400 c omplete eclipses 
seen from EXO 0748-676 (I Wolff et al.ll2009l '). the mid-eclipse 
times measured with RXTE-PCA are accurate enough to 
enable detection of very small orbital period glitches. 

Magnetic field cycling of the secondary star is as- 
sumed to be the likely cause for the observed orbital pe- 
riod glitches in EXO 0748-676 (|Wolff et al.ll2009l ). It is pro- 
posed that magnetic activity associated with the secondary 
star could be responsible for sudden changes in the orbital 
period. If the secondary star associated with XTE J1710- 
281 has strong, changing, magnetic activity, it can lead 
to changes in the structure of secondary star. A changing 
gravitational quadrapole moment can res ult into changes in 
the orbital period of the binary syst em (iLanza fc Rodong 
1 19991 : iTauris fc van den Heuvel|[200^ ). However, in case of 
XTE J1710-281, the o ptical counterpart has been discov- 
ered (|Ratti et al.ll2O10l ) but the type of the companion star 
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Figure 3. A sample of 2-20 keV folded light curves of XTE J1710-281. The light curves were folded with a period of 0.1367109674 d at 
an epoch of MJD 51250.924540. The normalised intensities during epoch-1 and epoch-2 have been rescaled, by adding constant numbers 
to the curves. The solid line in the middle light curve (epoch 2) shows the best fit five-parameter model to the X-ray eclipse. The best 
fit had a reduced of 689 for 650 d.o.f. The top panel shows the complete eclipse, whereas, the bottom panels show the ingress and 
egress of the eclipse phase. 



is not yet known. Therefore, it is difficult to make a state- 
ment on the probable cause for the changing orbital period 
of XTE J1710-281. 

We emphasise that if magnetic cycling of the binary 
components is indeed a reason behind the observed epochs 
of orbital period, then long term monitoring of XTE J1710- 
281, is required to determine the timescales of magnetic cy- 
cling of the secondary star. It may also be useful to foretell 
the distinct orbital period epochs of XTE J1710-281, if any. 
Lastly, the forthcoming Indian-satellite, ASTR OS AT with 
a very large area X-ray proportional counter (jPaull |2009| ) 
could be a boon in determining the orbital parameters of 
the system. 
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Table 1. Mid Eclipse Time measurements of XTE .11710-281. 



RXTE 


Orbital 


Mid eclipse time 


Icr Uncertai] 


Observation ID 


Cycle 


M,1D (d) 


(d) 


40407-01-03-00 


1 


51251.061141 


0.000006 


40135-01-40-00 


728 


51350.450003 


0.000035 


60049-01-01-03 


6452 


52132.983718 


0.000046 


60049-01-03-00 


7118 


52224.033183 


0.000012 


60049-01-04-00 


7567 


52285.416436 


0.000029 


60049-01-05-00 


7623 


52293.072234 


0.000006 


60049-01-06-00 


7667 


52299.087534 


0.000012 


60049-01-07-010 


9431 


52540.245711 


0.000064 


60049-01-07-01 


9432 


52540.382406 


0.000030 


80045-01-01-00 


11732 


52854.817626 


0.000006 


80045-01-01-01 


11733 


52854.954339 


0.000010 


80045-01-01-03 


11756 


52858.098681 


0.000009 


80045-01-01-04 


11757 


52858.235396 


0.000010 


80045-01-01-05 


11760 


52858.645512 


0.000009 


80045-01-02-00 


12432 


52950.515276 


0.000009 


80045-01-02-01 


12444 


52952.155841 


0.000012 


80045-01-02-03 


12445 


52952.292553 


0.000006 


80045-01-02-02 


12454 


529-53.522899 


0.000014 


80045-01-02-04 


12455 


52953.659644 


0.000019 


80045-01-02-05 


12456 


52953.796342 


0.000008 


80045-01-03-00 


12947 


53020.921420 


0.000012 


80045-01-03-01 


12957 


53022.288557 


0.000009 


80045-01-03-02 


12970 


53024.065797 


0.000017 


80045-01-03-03 


12971 


53024.202500 


0.000015 


80045-01-03-04 


12982 


53025.706300 


0.000012 


80045-01-03-05 


12983 


53025.843038 


0.000012 


80045-01-03-06 


12984 


53025.979734 


0.000008 


80045-01-03-07 


12985 


53026.116419 


0.000017 


80045-01-04-00 


15096 


53314.713316 


0.000030 


80045-01-04-01 


15107 


53316.217123 


0.000013 


80045-01-04-02 


15115 


53317.310808 


0.000017 


80045-01-04-05 


15129 


53319.224759 


0.000021 


80045-01-04-06 


15130 


53319.361463 


0.000021 


80045-01-04-07 


15131 


53319.498156 


0.000010 
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17718 


53673.169500 
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91045-01-01-01 


17719 


53673.306207 
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17740 


53676.177102 
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17741 


53676.313826 
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17785 
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0.000009 
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17786 


53682.465792 


0.000014 
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53683.969629 


0.000012 


91045-01-01-16 


17798 


53684.106327 


0.000010 


91045-01-01-17 


17799 


53684.243048 


0.000008 


91018-01-01-00 


20059 


53993.209860 


0.000021 
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20060 


53993.346545 


0.000012 
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54112.831940 
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91018-01-03-01 


21314 


54164.782093 


0.000006 
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22700 


54354.263475 


0.000022 


91018-01-07-04 


22758 


54362.192739 


0.000006 
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54407.033934 


0.000006 


91018-01-08-01 
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54407.170661 
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0.000012 


91018-01-08-03 
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0.000006 
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55208.980538 


0.000023 


94314-01-06-00 


30454 


55414.320400 


0.000030 


94314-01-06-00 


30455 


55414.457119 


0.000029 


94314-01-06-01 


30461 


55415.277413 


0.000017 



Table 2. Orbital epliemerides of XTE .11710-281. 



Pcirameter 


Best fit value from present analysis 


To (MJD) 


51250.924540 (4) 


Port (d) 


0.1367109674 (3) 


Porb (10 d d-1) 


-1.6 ^ Porb ^ 0.2 


Porb/Porb (10 * yr) 


-2.34 Porb/Porb ^ 18.7 
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